ensitivity, uncertainty and time dependence of arameters in a complex land surface model

This paper explores the use of Monte Carlo carbon cycle data assimilation within the generalized likelihood uncertainty estimation (GLUE) frameworktoevaluate thesensitivitiesofa wellknown complex land surface model (SiB v2.5) to its parameterization and the predictive uncertainty of simulated fluxes on a monthly basis, and for an entire year, at the WLEF talltower site in Park Falls, Wisconsin. An analysis is described wherein randomly generated parameter sets were ranked based on their capacity to simulate fluxes of latent (LE) and sensible heat (H) and the net ecosystem exchange of carbon (NEE) for each month of the year 1997 and for the entire year. Two criteria were used to evaluate the success of the simulations; the first evaluated the ability of SiB2.5 to simulate LE and H, the second included NEE as an additional constraint. The best-performing parameter sets for each criterion were used to assess model sensitivityto parameters, to calculate uncertaintybounds for predicted LE, H and NEE and to assess the information content of eddy covariance data on the analyzed time scales. Patterns in the sensitivity of the model to its parameterization and the uncertainty of the predictions were related to the physiological and phenological characteristics of the ecosystem, model structure and the relationship between deterministic models and comparatively stochastic measurements. The results show that model sensitivity varies through time for a larger set of parameters than those typically considered time varying in LSMs, thus optimization of model parameters on tower flux data should allow for variability at subannual time scales in order to capture the most information and best simulate fluxes. Further, constraining predictions annually versus monthly showed that some quantities (e.g. nighttime NEE) were on average better constrained annually, whereas other quantities that show more variability with vegetationphenology and structure (e.g.daytime NEEand LE) werebetter constrained monthly. The addition of the net ecosystem exchange of carbon in the data assimilation scheme improved model results by (a) constraining model parameterization (optimal parameter values), particularly during times of the year when the land surface was rapidly changing (spring and fall), and increasing the number of influential parameters, and (b) decreasing the uncertainty in NEE simulations (but not appreciably reducing uncertainty in LE andH simulations). Wealso foundthat therewasanirreducible levelofmismatchbetween the simulated and observed fluxes that could not be overcome through optimization due to variability in the observations and/or structural problems with the model. The uncertainty estimates can be used to characterize uncertainty in the simulations at multiple time scales. # 2007 Published by Elsevier B.V. avai lable at www.sc iencedi rec t .com journal homepage: www.e lsev ier .com/ locate /agr formet * Corresponding author. Tel.: +970-491-1964; fax: +970-491-1965. E-mail address: [email protected] (L. Prihodko). 168-1923/$ – see front matter # 2007 Published by Elsevier B.V. oi:10.1016/j.agrformet.2007.08.006 a g r i c u l t u r a l a n d f o r e s t m e t e o r o l o g y 1 4 8 ( 2 0 0 8 ) 2 6 8 – 2 8 7 269